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The  inertial  process  of  multiple  position  gyrocompassing  is  analyzed  to 
determine  performance  deviations  due  to  the  geokinetic  error  sources  of  tilt,  tilt 
rate  and  polar  wobble. 
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Fundamental  Geokinetic  Considerations 
in  Multiple  Position  Gyrocompassing 


1.  INTRODUCTION 


The  inertial  guidance  system  is  a  major  subsystem  in  many  Air  Force  weapon 
systems.  The  primary  purpose  of  the  guidance  system  is  to  provide  navigation 
and  other  pertinent  weapon  delivery  data  to  the  system  operators.  To  successfully 
perform  this  function,  it  is  necessary  that  the  inertial  guidance  system  be  accurately 
initialized.  Data  such  as  initial  position,  velocity  and  target  position  are  required 
as  inputs  to  the  navigation  computations.  In  addition,  the  inertial  reference  plat¬ 
form  must  be  aligned  with  respect  to  some  known  reference.  The  optical  alignment 
and  the  self-alignment  process  are  the  two  particular  techniques  most  frequently 
used: 

(1)  In  the  optical  alignment  technique,  the  inertial  platform  is  slaved  to  a 
known  optical  reference.  It  is  essential  to  have  a  known  stable  optical  reference, 
such  as  an  optical  cube  mounted  on  a  pier,  whose  azimuth  with  respect  to  celestial 
North  has  been  determined  through  extensive  measurement  to  an  accuracy  better 
than  the  azimuth  required  by  the  inertia’  reference  unit.  Additionally,  some  means 
must  exist  to  transfer  the  cube  azimuth  to  the  inertial  platform. 

(2)  In  the  self-alignment  process  technique,  the  inertial  platform  requires 
no  external  references.  Here  the  astronomic  azimuth  of  the  inertial  element  is 
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determined  by  making  measurements  of  the  Earth's  diurnal  rotation  vector,  W.  , 
and  the  Earth’s  gravity  vector  g. 

Under  Project  7628  (Geophysical  and  Geokinetic  Effects),  AFCRL  scientists 
have-  been  investigating  the  effects  of  geokinetic  disturbances  on  the  performance 
of  inertial  instruments  and  systems.  Since  inertial  instrumentation  requires 
electromechanical  components  whose  function  demands  high  sensitivity  to  mechani¬ 
cal  motions,  geokinetic  errors  sources,  if  not  properly  accounted  for,  will  obvious 
ly  degrade  performance  of  the  instrumentation.  The  purpose  of  this  report  is  to 
define  the  theoretical  performance  deviations  which  could  be  expected  in  a  self¬ 
alignment  process  for  inertial  instrumentation  called  "multiple  position  gyro- 
compassing.  "  These  deviations  are  a  direct  result  of  the  fundamental  geokinetic 
motions  of  tilt  and  polar  wobble. 


2.  MULTIPLE  POSITION  GYROCOMPASS  ANALYSIS 


In  multiple  positic.i  gyrocompassing,  a  precision  rate  integrating  gyroscope 
operating  in  the  rule  mode  is  sequenced  through  a  munber  of  different  attitudes 
relative  to  the  Earth  rate  vector  and  gravity  vector.  The  rate  measurements 
made  in  these  positions  are  then  differenced  and  averaged  to  obtain  a  bias-free 
astronomic  azimuth  of  the  sensors  reference  axis  that  is  independent  of  the  gyro 
drift  a>  d  bias. 

Figure  1  is  a  pictorial  representation  of  a  single  degree-of -freedom  rate 
integrating  gyroscope.  When  operating  in  the  rate  mode,  the  gyro  performance 
model  is  given  by: 


W  .  WInA  +  B  !m  hPIM  fSRA 


where 

WIRA  =  inPut  refe}  ence  axis  component  of  the  angular  velocity  of  the 
gyro  case  with  respect  to  inertial  space. 

B  =  the  constant  portion  of  the  gyro  output,  bias. 
a  =  a  random  uncertainty  of  the  gyro  output. 

PgRA  3  pendulosity  (mass  unbalance  times  distance)  along  the  spin 
reference  axis. 

PjrA  =  pendulosity  along  the  input  reference  axis, 
fjRA  =  specific  force  along  the  input  reference  axis. 
fgjjA  “  specific  force  along  the  spin  reference  axis. 

H  =  angular  momentum  of  the  gyro  rotor. 
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Figure  1.  Single  Degree-oi -Freedom  Rate  Integrating  Gyroscope 


If  such  a  gyro  is  sequenced  through  four  positions  given  below,  it  is  possible  to 
determine  the  astronomic  azimuth  of  the  sensors  spin  reference  axis  from  the  rates 
measured  by  the  gyro. 

Position  1 

IRA  East 

ORA  Down 

SRA  North 

where 

Down  =  along  the  direction  g. 

North  =  in  ?  plane  perpendicular  to  g  in  the  direction  of  the  local  horizontal 
projection  of  the  principle  axis  of  the  Earth's  figure. 

East  s  forme  a  right  handed  orthogonal  set. 

The  nominal  position  of  the  gyro  frame  (g')  is  defined  by  the  North,  East,  and 
Down  reference  directions  as  defined  above.  The  actual  position  of  the  gyro  frame 
(g)  will  differ  from  the  nominal  gyro  frame,  because  the  gimbal  system  which 
supports  the  gyroscope  will  not  be  perfectly  aligned  initially  and  will  continuously 
change  orientation  because  of  clianges  in  the  gimbal  base  with  respect  to  the  gravity 
vector  or  because  of  changing  environmental  conditions.  In  either  case,  the 


Position  2  Position  3  Position  4 

IRA  West  IRA  East  IRA  West 

ORA  Up  ORA  Up  ORA  Down 

SRA  North  SRA  South  SRA  South 


attitude  of  the  g  frame  with  respect  to  the  g'  frame  can  be  determined  from  th- 
direction  cosine  matrix 


where  the  ti’e  angles  e^,  eE,  and  eD  have  been  considered  to  be  small  enough  ao 
that 


cos  e.  *>  1 

J 

ainej  *  cj 


and 


wg, 

g  g 


With  these  definitions  and  the  assumption  that  the  gyrocompass  is  stationary  with 
respect  to  a  point  on  the  surface  of  the  Earth,  it  is  possible  to  coordinate  the 
specific  force  vector  and  the  Earth's  daily  rotation  vector  in  the  actual  gyro  frame. 

In  the  case  where  nominal  motion  of  the  gyroscope  with  respect  to  the  mass 
center  of  the  earth  is  zero,  the  specific  force  vector  and  the  gravity  vector  are 
identical.  Then  the  specific  force  can  be  written  as 


The  Earth's  angular  velocity  can  be  coordvnatized  in  the  Earth  centered  inertial 
frame  as 


where 

is  the  primary  component  of  Wj  ,  and 
v»x,  W  are  equatorial  components  ot  Wjff. 


^^v?^?*<?^?.*>''-’C*r^>*>*^r»'-w’r.' . . .  . . .  .  „  ,  . 


tt  *  ’tj;  ,  vtt  - 


Then  the  Earth's  angular  velocity  in  the  g'  frame  is  given  by 


wf'  = 

-ie 


-  sin  L  cos  X  V,JX  -  sin  L  sin  X  Wy  +  cos  L  W„ 

-  sin  X  V/  +  cos  X  W 

x  y 

-  cos  L  cos  X  Wx  -  cos  L  sin  X  Wy  -  sin  L  W„ 


where 

L  is  the  reduced  astronomic  latitude 
X  is  the  reduced  celestial  longitude 
X  =  t  +  wzt 

l  —  reduced  terrestrial  longitude  from  Greenwich 
t  —  is  time 

The  earth's  angular  velocity  coordinatized  in  g  frame  then  becomes 


cos  L  W  +  sin  LW  -  sin  L  cos  XWv  •  sin  L  sin  X  w 
z  ^  l  x  y 

W’  =  {  -  sin  X  +  cos  X  Wy  -  eD  ccs  L  W2  -  e N  sin  L  Wz 

-sin  L  W  +  cos  L  W  -  cos  L  cos  X  w„  -  cos  L  sin  X  w 
z  z  iii  x  y 


’ie 


From  the  fundamental  law  of  vector  addition  of  angular  velocities,  the  gyro  case 
angular  rate  can  be  written  as 

Wg  *=Wg  +  wg  +  w*\ 

-Ig  -ie  -eg'  -g'g 

or 


W. 


•ig 


cos  L  Wz  +  sin  L  WzCg-  sin  Lcos  X  Wx  -  sin  L  sin  XWy  +  eN 

-sin  X  Wx  >  cos  X  W^  -  cos  L  e^  -  sin  LW^^  +  eE 

-3in  LW  +  cos  LW  e.-  -  cos  L  cos  X\V  -  cos  L  sin  XW  +  6^ 
z  z  £i  x  y  l/ 


(4) 


With  the  relationships  Eqs.  (4)  and  (3)  substituted  into  Eq.  (1)  for  the  positions 
stated  in  Eq,  (2),  the  following  four  indicated  rates  are  obtained: 

W,  *  -sin  XWv  +  cos  X  W  -  e  cos  LW„  -  eN  sin  LW  +  B  +o. 

1  x  y  D  z  «  z  J 

+  n  (PSRA  €N  g  +  PIRA  CE  g)  +  f  E 
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W2  =  +sin  X  Wx  -  cos  X  Wy  +  eD  cos  L  Wz  +  sin  L  W2  +  B  +  Og 

+  H  (_PSRA  eNg  +  PIRA€Eg>  ‘  eE 
Wg  *  -sin  X  Wx  +  cos  X  Wy  -  cos  LWz  -  sin  L  Wz  +  B  +  Og 

+  R  (+PSRA  eNg  ‘  PIRA  eEg)  +  CE 
W4  =  +sin  X  Wx  -  cos  X  cos  L  Wz  +  sin  L  W2  +  B  + 

+  H  ("PSRA  eNg  ‘  PIRA  CEg)  ‘  eE  ' 

The  quantity 

W  =  1/4  (Wj  -  W2  +  Wg  -  W4) 
is  formed  and  is  seen  to  be  equal  to 

W  =  -eD  cos  L  Wz  -  e^j  sin  L  V>’2  -  sin  X  +  cos  X  \V^ 

+  PSRA  cNg  +  °  +  £E 

v  here 


In  addition  it  should  be  noted  that  ar  y  mass  unbalance  drift  terms,  due  to  the 
specific  force  component  along  the  output  axis,  will  also  be  eliminated  by  the 
averaging  process  above. 

The  normal  calculation  made  is 

-W 

£D  cos  L  VVz 

V’hic'n  is  the  astronomic  azimuth  of  the  gyro  spin  reference  axis  and  is  in  error  by 
the  four  following  quantities: 


sec  L  L 
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sec  L  a 
2  “  Wz 

/  ,  „  PSRA g  \  _ 

e3  =  ^  tan  L  -  sec  L  jcN 

.  S  W  X  W 

sin  a  x  _  cos  A  y 

4  cos  L  W_  cos  L  W 
z  z 

which  represent  the  errors  due  to  the  imperfect  level  of  the  instrument  about  the 
North-South  axis,  th-  tilt  rate  about  the  West-East  axis,  the  component  of  the 
Earth's  rotation  vector  in  the  reduced  equatorial  plane,  and  the  statistical  behavior 
of  the  gyro  drift. 

3.  PERFORMANCE  DEVIATIONS  Dl.'E  TO  GEOKIN ETIC  EFFECTS 

Of  the  four  error  sources  ej  through  e4  defined  in  Section  2,  three  are  due  to 
geokinetic  disturbances.  Before  further  investigation  of  these  deviations,  however, 
we  will  comment  on  the  effect  of  the  random  fluctuations  of  gyro  drift. 

In  general,  the  statistical  description  of  gyro  drift  is  not  agreed  upon. 
Weinstock  (1964)  has  shown  large  variations  in  the  statistical  descriptions  of 
identically  designed,  high-precision  inertial  gyros.  The  model  adaptea  here  is  in 
agreement  with  the  majority  of  analytical  attempts  to  describe  the  gyro  drift.  The 
drift  model  to  be  used  is 

o'  =  -  l/k  o'  +  n(t) 

where 

o’  =  the  random  gyro  drift  rate, 
k  =  the  correlation  time  of  the  random  process. 
n(t)  =  gaussian  distributed  white  noise. 

The  gyro  drift  uncertainty  is  modeled  as  a  Gaussian  Markov  process  exponentially 
correlated  with  coi  relation  time  k.  The  mean  of  the  gyro  drift  at  any  time  t  is 
given  by 

<  o'(t)>  =  e  k  <o'(to)> 


kV'WiJjwway*  „ _ _ 
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and  th<  variance  is  given  by 

/  !±V 

E0,fO  ■  ^  ^1  -  e  k 

where  N  is  the  constant  value  of  the  power  spectral  density  of  the  noise,  n(t). 
For  time  periods  (t-tQ)  small  compared  to  the  correlation  time,  the  variance  of 
the  gyro  drift  is  given  by 

Effi(t)  =  (t-t0)N 

or  the  process  is  similar  to  a  random  walk.  For  time  periods  longer  than  the 
correlation  time,  the  gyro  drift  is  bounded  in  that 


Britting  et  al  (1971)  gives  the  following  values  for  N  as  a  function  of  when  xhe 
instrument  was  designed. 

N  in  min2 /hr3 

Production:  0.0735 

State  of  the  Art:  0.0368 

Future:  0.00735 

The  random  fluctuation  (a)  which  exists  at  the  end  of  the  differencing  and  averaging 
process  can  be  determined  as  a  function  of  the  gyro  statistical  description.  The 
mathematical  expectation  of  a  is  given  by  ' 

f 

<  a>  =  1/4  <  -  a2  +  a„  -  a4  >  . 

If  the  time  interval  (t^-tg)  is  small  compared  to  the  correlation  time  k,  then 
<al>  *  <a2>  *  <a3>  *  <  V 


and 


<  a>  =  0  . 


This  condition  is  realistic  since  in  a  typical  multiple  position  gyrocompass 
mechanization,  the  measurement  interval  (t^  -tQ)  is  approximately  10  minutes. 


while  tiie  gyro  correlation  time  is  about  24  hours  (Blitting  et  al,  1971).  The 
variance  of  a  is  given  by 


Ea=<a2>  Iii<(0l'02  +  a3-a4)2> 


Eo  =  lV  <  al  +  a22  +  °3  +  °42  > 


+  8  <  °1°3  +  ff3a4  *  ala2  -  °1°4  -  a2a3  ’  a3a4  >  * 

Since  the  measurement  period  is  small,  each  value  of  the  drift  rate  will  be  indepen¬ 
dent  of  the  next  value.  With  additional  manipulation  the  variance  of  a  becomes 


E  *  rE„  +  E  +  E„  +  E„  1 
a  16  [  (Jj  o2  o3  cta  j 


Then,  using  the  expression  for  the  variance  valid  for  small  time  intervals 

■  -nr N  W  • 

Using  the  stated  values  above  of  N,  the  variance  of  the  output  of  the  multiple  posi¬ 
tion  gyrocompass  indication  at  45°  latitude  is  given  below. 


Class  of  Gyro 
Production : 
State  of  the  Art: 
Future: 


9.  0  sec^ 
4.  5  sec2 
0.  9  sec2 


In  all  cases  it  can  be  seen  that  the  random  fluctuations  of  the  gyro  drift  are  a  zero 
mean  variable  but  can  contribute  significantly  to  the  variance  of  the  system  output 
with  the  best  case  (future  gyro,  3-standard  deviations)  having  a  value  of  2.  9  sec. 

The  error  due  to  the  tilt  about  the  North-South  direction,  as  defimd  above,  has 
two  contributors  one  of  which  i|  due  to  the  imperfect  performance  of  the  gyro.  Note 
that  this  error  could  be  tuned  out,  at  least  in  theory,  by  fabricating  a  rate  integrat¬ 
ing  gyro  with  the  spin  mass  unbalance  satisfying  the  following  relationship: 

W,  sin  L  H 

T".  _  Z _ 


This,  however,  is  impractical  for  several  reasons.  The  contribution  due  to  this 
error  source  is  small  in  that  for  a  typical  large  mass  unbalance,  the  error  is  on 
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the  order  of  10  meru/g*.  Hence  at  45°  latitude 
-  0.014 

which  is  second  order  at  the  chosen  latitude.  Hence  the  major  error  source 
associated  with  North-South  tilt  is  due  to  the  primary  geometric  source,  tan  L. 

At  45°  latitude,  every  arc  sec  of  tilt  results  in  a  deviation  of  one  arc  sec  in  the 
gyrocompass  output.  Note  that  tilt,  as  defined  here,  is  a  rotation  of  the  case  of 
the  gyrocompass  with  respect  to  g,  the  gravity  vector. 

The  deviation  in  the  gyrocompass  output  due  to  the  tilt  rate  about  the  East- 
West  line  has  a  large  sensitivity  in  that,  at  45°  latitude, 

ej  =  292  sec 

where  the  East  tilt  rate  is  given  in  meru.  The  largest  tilt  rates  recorded  at 
AFCRL's  Haskell  Observatory  were  4  sec/day,  or  3.  6  X  10-3  meru  (Cabaniss, 
1972).  The  corresponding  deviation  in  gyrocompass  output  would  be 

*  1.  08  sec 

which  again  is  small.  However,  since  the  sensitivity  of  the  output  to  East  axis  tilt 
rate  is  relatively  large,  provision  should  be  made  to  monitor  the  East  axis  tilt 
rate  and  compensate  the  output  if  required. 

The  remaining  error  term  is  due  to  ihe  physical  nature  of  the  rotation  of  the 
Earth.  With  a  simple  model  of  the  Earth's  polar  motion,  the  equatorial  compon¬ 
ents  of  the  Earth  rotation  vector  \an  1  -  written  as 

W  *  cos  (Wot  +  4E>) 

X  6  bl 

W  *  W'e  sin  (Wgt  +  <t> ) 
where 

Wg  —  magnitude  of  the  equatorial  component  of  W.g. 

Wg  —  angular  rate  of  rotation  of  the  equatorial  components  about  the 
principal  Z  axis  of  the  Earth. 

O  —  is  phase  angle  relating  position  of  instantaneous  pole  at  t  =  tQ. 


*  1  meru  =  earth  rate  -  0.  015°/hr 
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With  this  representation  —  which  assumes  mean  pole  position  of  the  principal  axis 
of  the  Earth  —  the  gyrocompass  error  due  to  polar  motion  can  be  expressed  as 


W. 


cos  LW_ 


cos  (Wgt  +  <E>  3 ) 


or 

*  -  a  sec  L  cos  (Wgt  + 

where  4>j  differs  from  3>  by  the  terrestrial  longitude,  and  a  is  the  distance  of  the 
instantaneous  pole  from  the  mean  pole  expressed  in  terms  of  central  angle 
(1  sec~  100  ft).  At  45°  latitude  the  maximum  deviation  observed  due  to  polar 
motion  effects  would  be 

eg  (max)  =  a  yr. 

This  functional  dependency  has  been  derived  for  a  multiple  position  gyrocompass, 
but  can  be  shown  to  be  valid  for  all  other  types  of  inertial  system  gyrocompassing. 

Additional  geodetic  errors  may  exist,  but  are  only  due  to  incorrect  input  data, 
analysis,  or  operation,  rather  than  the  correct  definition  of  the  physical  vectors 
being  measured.  For  instance,  if  the  geographic  latitude  is  used  in  the  calculation 
of  the  azimuth  data,  rather  than  the  reduced  astronomic  latitude,  an  additional 
error  will  exist  equal  to 

e5  =  CD  ton  L  6 

where  6  is  the  difference  between  the  geographic  latitude  and  the  reduced  astro- 
nomic  latitude  expressed  in  radians.  This  again  is  an  extremely  small  effect.  The 
major  geokinetic  disturbance,  which  has  not  been  discussed  here,  is  seismic, 
motions.  In  general,  these  motions  will  have  an  effect  on  the  indicated  azimuth  in 
that  these  motions  degrade  tne  gyro  such  that  the  gyro  performance  model  given 
above  is  not  completely  applicable.  The  explicit  relationships  between  seismic 
motions  and  gyroscope  performance  is  the  subject  of  continuing  AFCRL  research 
and  will  be  reported  on  in  the  future. 


4.  SUMMARY 

It  has  been  shown  that  the  inertial  process  of  multiple  position  gyrocompassing 
is  subject  to  deviations  attributable  to  the  geokinetic  motion  of  North  axis  tilt. 


East  axis  tilt  rate,  and  polar  motion.  The  process  has  highest  sensitivity  to  East 
axis  tilt  rate.  The  tilt  about  the  North  axis  effects  the  output  directly  with  a 
sensitivity  of  tan  L.  Deviations  in  performance  due  to  polar  motion  are  small, 
but  probably  within  the  measurement  capability  of  today's  instrumentation. 
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